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ABSTRACT 

The natural approach to a spectral analysis of data distributed on the sky employs 

spherical harmonic decomposition. A common problem encountered in practical astronomy 

ON is the lack of full sky coverage in the available data. For example, the removal of the galactic 

G\ 

plane data from the COBE DMR sky maps compromises Fourier analysis of the cosmic 

H 

microwave background (CMB) temperature distribution due to the loss of orthogonality 
of the spherical harmonics. An explicit method for constructing orthonormal functions on 
an incomplete (e.g. Galaxy-cut) sphere is presented. These functions should be used in 
the proper Fourier analysis of the COBE DMR sky maps to provide the correct input 
O for the determination of the spectrum of primordial inhomogeneity. The results of such 

an analysis are presented in an accompanying Letter. A similar algebraic construction of 
appropriate functions can be devised for other astronomical applications. 
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1. INTRODUCTION 

Analysis of the first year of COBE DMR mission data resulted in the discovery of 
temperature anisotropy on the microwave sky of cosmological origin (Smoot et al. 1992, 
Bennett et al. 1992, Wright et al. 1992). A summary of the results from two years 
of observations was presented in Bennett et al. (1994). The COBE DMR instrument 
has observed the microwave sky at three frequencies (31.5, 53, and 90 GHz) with pairs 
of 7° (FWHM) antennas separated by 60°. The temperature difference measurements 
were used to construct whole sky anisotropy maps (Janssen &: Gulkis 1992) binned into 
6144 equal area (~ 2°6 X 2°6) pixels on a quadrilateralized spherical cube projection 
(White &: Stemwedel 1992). The temperature fluctuation assigned to a given pixel p, 
A(p) = ST /T(8 p ,(j)p), is primarily due to cosmic microwave background (CMB) radiation, 
emission from the Galaxy, and receiver noise. Systematic errors in the measurements have 
been carefully studied (Kogut et al. 1992, Bennett et al. 1994), and are ignored in the 
present analysis. 

The observed full sky temperature distribution, A(p) = Acmb(p) + ^Galaxy(p) + 
^■Noise{jp)i can be expanded in spherical harmonics 

*-m ax 

a/„\ [\CMB i Galaxy] DMR , Noise \ „„pix v in 1 \ 

A (P)=Z^ 1^ \[ aim + \ W t +a t™ ] W t Y £m (0 p ,(l) p ). (1) 

1=0 m=-l 

Here, the Yi m are real, orthonormal spherical harmonicsf, wf MR are the filter coefficients 
of the DMR beam pattern (Wright et al. 1994), and w^ iX are the filter coefficients used 
to model the extra smoothing due to pixelization of the data (actually a circular top- 
hat window of solid angle Cl p i x = 47r/6144). ai m are harmonic coefficients of the CMB 
anisotropy, galactic emission, and instrument noise. In the sum over £, the lower limit is 
taken deliberately at £ = 0; although the instrument does not directly measure the absolute 



f Y tm {6A) = V^VfeS^™'^ 08 ^^^' Where m = (( 2?r )" 1/2 ' ^ 1/2 cosm0, 
or tt" 1 / 2 sin \m\<f>), for m = 0, > 0, or < 0. P™ (x) m =° (-l) m (I - x 2 ) m / 2 P^x), where 
Pt(x) = ^ j-[{x 2 - l) 1 are Legendre polynomials. 
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mean temperature on the sky, both a 00 a axy and an arbitrary (but small) offset introduced 
by the map making algorithm exist. Similarly, the a priori CMB dipole anisotropy removal 
employed in the two year data analysis (Bennett et al. 1994) is only approximate and leaves 
a small residual 1 = 1 component in the maps. Due to the pixel size, a formal upper limit 
in the sum over I is near I 

max ~ 100 for the raw, noise-dominated maps. However, the 7° 
FWHM DMR antennas are sensitive to sky signals only up to I ~ 30, so that at I ^ 30 
A(p) in eq.(l) is, in practice, dominated entirely by the a^ r ^ lse term. 

A major aim of the DMR data analysis is to extract the cosmological signal A.cmb(p), 
or the coefficients afJf B from the sky maps. To achieve this, one needs to identify and 
separate the signals due to receiver noise and galactic emission. The receiver noise in 
a given map can be well characterized statistically using the known properties of the 
radiometer and the sky coverage. The dominant galactic emission can be eliminated most 
easily by removing all pixels within 20 degrees of the galactic plane. A more detailed 
approach is to model the galactic emission in detail utilizing the multifrequency nature 
of the data (Bennett et al. 1992). The results of such modelling indicate that galactic 
emission at high latitudes is relatively weak, while the emission in the plane cannot be 
removed accurately enough to warrant using the data for anisotropy studies. Thus the 
technique of cutting the plane and ignoring the high-latitude emission appears to yield 
the best combination of minimizing galactic emission while maximizing the signal to noise 
ratio, since the modelling technique includes the relatively noisy 31.5 GHz DMR data. 

Implementing a galactic cut in the data creates a problem for Fourier analysis of the 
sky maps, since the drastically reduced sky coverage (e.g., a |fe| = 20° Galaxy cut results 
in ~ 35% loss in the solid angle) does not provide support for the orthogonality of the 
spherical harmonics (eq. 1). 

In this Letter, a method for generating an orthonormal basis of functions on the pix- 
elized, Galaxy-cut COBE DMR sky is presented. This method is a practical implemen- 
tation of the Gram-Schmidt orthogonalization (see e.g. Kostrikin &: Manin 1989) of a set 
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of linearly independent vectors. These algebraically constructed orthonormal basis func- 
tions are used for (1) a proper Fourier analysis of the reduced coverage sky maps, (2) a 
derivation of the correlation matrices of Fourier amplitudes for the model anisotropy and 
the instrumental noise pattern, and (3) a computation of the likelihood function for the 
cosmological structure formation model parameters to be inferred from the experimental 
data. 

Hereafter: bold, upper case letters denote matrices; bold, lower case letters denote vec- 
tors; subscripts i, j, k, I denote indices of orthonormal functions; and indices or subscripts 
p, q are pixel labels, which identify the angular position on the sky. 

2. ORTHONORMAL FUNCTIONS AND FOURIER ANALYSIS 
ON AN INCOMPLETE SPHERE 
Consider a linear space Ti max of bounded functions on a sphere spanned by spherical 
harmonics up to order l max . Let us form an (! moi + l) 2 -dimensional vector y = (Yo,Oj 
Yi-i, Yi, , Yi,i, Y 2 - 2 , Y 2 -i, Y tmax ,t max ). An index pair, (£,m), of the spherical 
harmonic degree and order is mapped into i = £ 2 -\-£-\-l-\- m, a single index of the vector 
yi. The inverse mapping is given by I = integer — 1), and m = i — (£ 2 + £ + 1). A 
function / £ J~i max is defined as a finite Fourier series = ^(^ aa: + 1 ) a^y^O) = a T -y, 

where a,- = J" 47r dft, yj-(O). For any functions /, g £ J~i max define the scalar product 
specific to the pixelized sky: 

(f9){f.Usk 9 }=^P» E M M ™ J & W) g{Q), (2) 

pE{full sky} 

where {full sky} denotes the entire set of 6144 sky map pixels. Spherical harmonics are 
orthonormal in the actual DMR pixelized sky scalar product, • y T ) ^ ull sky ^ = I, to an 
absolute accuracy better than 10 -3 for £ < 30. 

When a Galaxy cut is applied to the sky map the set of pixels that enter the scalar 
product is reduced. An explicit notation for this is 

U 9) {cut sk y} = Yl fiP) 9(P)- (3) 

pE{cut sky} 
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For example, the Galaxy cut at |fe| = 20° reduces the number of pixels from 6144 to 4016. 
If only a subset of a sphere provides support for the scalar product, then the spherical 
harmonics are linearly independent due to their different functional forms, but are not 
orthogonal. A new functional basis, orthonormal in the scalar product of eq. (3), can be 
algebraically constructed from the functions y as follows. 

Evaluate the coupling matrix for the spherical harmonics on the cut sky 

(y - y ) {cut sk v } = ( 4 ) 



pix 



Here, y are spherical harmonics with amplitudes slightly suppressed by the factor w ^ 
(e.g. w^ 1 q ~ 0.93), due to pixelization. The Kowalewski-Gram determinant det W (San- 
sone 1959), is an indicator of whether a set of underlying functions is linearly independent. 
In the case of reduced sky coverage, W is a non-diagonal, symmetric, non-singular, pos- 
itive definite matrix. Hence, W can be Choleski-decomposed into the product of a lower 
triangular matrix L and its transpose 

W = L L T . (5) 

The essential point of this method relates to the observation that Choleski decomposition 
of the coupling matrix W provides a direct means of Gram-Schmidt orthogonalization. 
Compute the inverse matrix Y = L _1 , and define a new set of functions on the cut sky 

*(p) = r-y(p). (6) 

By construction 

, .* T \ =r ./ y .yT\ . r T = r • l • l t • r T = i. (7) 

Y /{cut sky] W 3 /{cut sky} V ; 

Hence, the functions ij) are orthonormal on the pixelized, cut sky, and are suitable for the 
Fourier analysis of functions in Ti max with support restricted to the cut sky. This algorithm 
has been numerically implemented for l moi = 30. Orthonormality of the resulting 961 ij) 
functions was achieved to an absolute accuracy better than 10 -7 . 



Using the new basis ij) , a function f{jp) = a T • y(p) can be Fourier decomposed in the 
cut sky: 

f(p)\ P e{cut sky} = c T i> (p), where c = (/ ^ ) {cut sky} . (8) 

It follows from eq.(6) that the relations between expansion coefficients in the y and ij) 
bases are 

c = L T • a, and a = T T ■ c. (9) 

It is important to note that due to the triangular form of L and T, each coefficient c ? - is a 
combination of coefficients aj of order j > i. Hence, truncating the first four components 
of the vector c unambiguously removes the monopole and dipole degrees of freedom from 
the temperature map, but does not affect the information content of the map for higher 
orders of I. This property is a consequence of the specific ordering of the y basis. 

3. MODEL CMB ANISOTROPY CORRELATIONS IN THE i/> BASIS 
Mainstream cosmological structure formation models invoke homogeneous, isotropic, 
Gaussian fields of small amplitude curvature perturbations superposed on a flat back- 
ground, specified by the spatial correlation function, or its Fourier transform, the power 
spectrum. The angular distribution of CMB temperature anisotropy induced by such cur- 
vature perturbations (Sachs &: Wolfe 1967) is easily described in Fourier language (e.g. 
Peebles 1981). Individual spherical harmonic coefficients of such a CMB anisotropy field 
are Gaussian-distributed in the theoretical ensemble of initial conditions for structure for- 
mation. The variances of the probability distributions of individual modes, ™)}' 
are uniquely expressed as integrals over the power spectrum, and depend only on £ due 
to the statistical isotropy of the CMB temperature field. Moreover, the ensemble corre- 
lation matrix of the coefficients slcmb is diagonal: (slcmb • & cmb) = diag |^ a f(^ ™)}}" 
Using eq.(9) one readily obtains the correlation matrix for the Fourier coefficients of the 
theoretical cut sky CMB anisotropy 

CcMB = (ccMB ■ Ccmb) = ^ T " (&-CMB " &CMb) " ^; (10) 
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where the effect of DMR beam smearing is included in the vector slqmb with elements 
-CMB _ a CMB W DMR ^ an( j pi xe li za tion smoothing is included in the matrix L (eqs. 4, 
5). The Gaussian ensemble probability density of cqmb is specified by the matrix Cqmb- 
P(c C mb) dccMB oc dccMB exp (-c£ MB • C'^ 1 MB ■ c C mb/2) / \/det(C cms) . 

In the long wavelength perturbation regime probed by COBE DMR, the power spec- 
trum of density perturbations is usually assumed to be a power law, P(k) oc k n . From this 
and previous assumptions, one can derive a compact formula for the expectation values 
of theoretically predicted CMB anisotropy multipole amplitudes (see Bond &: Efstathiou 
1987, and Fabbri, Lucchin, &: Matarrese 1987) 

\ 2 r (/ + n=i) r (j^) 

a i{l,m)j - a 2 Y(l+5=™') T ' ^ ' 

where a 2 is the rms quadrupole coefficient, often specified as Qrms—PS — "^^2 ■> and 
r is the gamma function, not to be confused with the matrix T. A broad category of 
inflationary structure formation models predicts n ~ 1, the Harrison- Zel'dovich spectrum. 
Using existing observational data to test the hypothesis that eq.(ll) provides a satisfactory 
description of the universe on large scales is a critically important aspect of contemporary 
cosmology. 

4. NOISE CORRELATIONS IN THE i/> BASIS 
The effect of radiometer noise on the sky maps constructed from the COBE DMR 
observations can be described to leading order (i.e. neglecting the small noise correlations 
between pixels separated by 60°, see Lineweaver et al. 1994) as a spatially uncorrelated 
Gaussian process of vanishing mean, {A N (p)) ^ noise ensem u e } = 0' an ^ P lxe ^ dependent 
variance (A N (p)A N (q)) {noise ensemble} = <T 2 (p)6 pq . The function <r 2 (p) is determined from 
the known noise characteristics of each individual radiometer channel and the number of 
observations per pixel. cr(p) varies with both frequency and channel, and combinations 
thereof, and can be determined to a few percent accuracy. 
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Fourier decomposition of the random noise pattern on the cut sky yields 



^N(p)\ P e{cut sky} = Cn ■1> (p), where c N = (A N i/> ) {cut sky} . (12) 

On averaging over an ensemble of Gaussian noise realizations, (cjv) , n ■ c ori oc . , x = 0, 
and the correlation matrix of the Fourier components of the noise pattern is 



C;v = (cjv 



{noise ensemble} 

(A N (p) &N(q)) {noise ensemble} 4 (p) $ T (q)) ( U ) 



= Vpix (<T 2 (P) 4 (P) ^ T (P)] 

\ ' pE.{cut sky} 

Analogously, it can be shown that 

(ciCjC k ci) {noiae ensemble} = (ciCj) (c k ci) + (c iCk ) (cjCl) + (aci) (cjc k ) , (14) 

and similarly for higher even moments, and that all odd moments are identically zero. 
Hence the probability distribution of the noise multipoles cjv is purely Gaussian, fully 
specified by the matrix Cjy: P(cjv)c?cjv ex dc^ exp (— • C^ 1 • cjv/2) / «J det(Cjy) . 

Spatial uniformity in the noise distribution, cr(p) = const, would result in a diagonal 
correlation matrix Cjy. The actual non-uniform noise distribution in the COBE DMR 
sky maps is specified by the strongly diagonally dominated, symmetric matrices Cjy. 

5. LIKELIHOOD ANALYSIS 
Since Gaussian statistical ensembles were invoked in §3 and §4 for the modelling of 
both cosmological and receiver noise temperature fluctuations in the cut sky maps, the 
exact probability distribution for the measurable Fourier amplitudes, c = cqmb + cjv, of 
the noise contaminated cosmological CMB anisotropy on the cut sky is 

J2. p — |c T -(C C Mfl+Cjv) _1 -c 

P(c)dc = , , r/ , , (15) 

The monopole and dipole components of the sky maps, which are not physically relevant 
to the power spectrum estimation, can be exactly removed from the analysis by integrating 
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over the first four components of the vector c. This operation removes the first four rows 
and columns of the correlation matrix, and reduces its rank to N = (£ max + l) 2 — 4. The 
remaining N elements of the vector c retain the entire information content of the sky map 
within the spectral range £ £ [2, £ m ax]- The need to choose a specific value for l moi in a 
practical implementation of a Fourier decomposition of the sky maps creates the problem 
of leakage of signal from the spectral range £ > l moi into the £ £ [2,£ max \ regime. Clearly, 
it is impossible in the analysis of an unfiltered sky map at a single frequency to separate 
correctly the true sky signal from the high-^ noise contamination. However, the definition 
of the probability distribution of the measurable Fourier amplitudes c in eq. (15) can be 
naturally extended to a multifrequency analysis. For example, a simultaneous analysis of 
the DMR 53 GHz and 90 GHz data involves a composite correlation matrix of the form 

^ _ / CcMB + C./V53 CcMB \ ,-, a s 

^53090 — I n n , n I , (,-LDj 

and a measured amplitude vector 

^53090 = (^53) ^9o) ■ (17) 

The cross-correlation terms in eq.(16) limit the influence of the high-^ frequency dependent 
noise leakage (and, in fact, of any frequency dependent signal contamination due to the 
galactic foreground emission or systematic effects) on the inferred frequency independent 
cosmological signal. In a simultaneous analysis of multifrequency data sets full advantage 
is taken of both the auto-correlation structure of individual maps and the cross-correlation 
between different frequency maps. Thus, all of the available information from the multiple 
experimental data sets is utilized. 

By assuming values for the cosmological parameters a,2 and n (eq. 11), one can use 
eq.(15) in a frequentist Monte Carlo simulation of synthetic vectors c. However, since we 
know the measured vector c and want to know the likely values of the ensemble averaged 
power spectrum parameters CL2 and n, which specify the matrix Ccmb-, a Bayesian ap- 
proach (e.g. Berger 1980) is suggested, wherein eq.(15) is a definition of the likelihood 
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for the parameters a,2 and n, given the specific data set. This approach and subsequent 
analysis is discussed in an accompanying Letter (Gorski et al. 1994) on the application of 
the formalism presented in this work. 

6. SUMMARY 

An algebraic method for constructing orthonormal basis functions on an incomplete 
sphere, designed for use in the Fourier analysis of the Galaxy-cut COBE DMR sky map- 
s, has been presented. The method involves a completely general implementation of the 
Gram-Schmidt orthogonalization of linearly independent vectors, and could be easily mod- 
ified for use in other astronomical projects with incomplete sky coverage. 

The primordial power spectrum determination based on the likelihood method present- 
ed in §5 applied to the two year COBE DMR sky maps is presented in an accompanying 
Letter (Gorski et al. 1994). 

I am thankful to A.J. Banday, C. Bennett, G. Hinshaw, A. Kogut, and E.L. Wright for 
numerous discussions during the development of this project and help in improving the 
manuscript. 
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